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The thermal equilibrium of string gas is necessary to activate the Brandenberger-Vafa mechanism, 
which makes our observed 4-dimensional universe enlarge. Nevertheless, the thermal equilibrium is 
not realized in the original setup, a problem that remains as a critical defect. We study thermal 
equilibrium in the Hagedorn universe, and explore possibilities for avoiding the issue aforementioned 
flaw. We employ a minimal modification of the original setup, introducing a dilaton potential. Two 
types of potential are investigated: exponential and double-well potentials. For the first type, the 
basic evolutions of universe and dilaton are such that both the radius of the universe and the dilaton 
\^ • asymptotically grow in over a short time, or that the radius converges to a constant value while the 

I dilaton rolls down toward the weak coupling limit. For the second type, in addition to the above 

, solutions, there is another solution in which the dilaton is stabilized at a minimum of potential and 

' the radius grows in proportion to t. Thermal equilibrium is realized for both cases during the initial 

. ' phase. These simple setups provide possible resolutions of the difficulty. 

I ^ . PACS numbers: 98.80.Cq, 11.25.Wx 
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String theory represents the most useful and promising candidate for a unified theory of the fundamental interactions 

■ including gravity. From cosmological aspects, an ultimate goal of such theory is to resolve three interesting problems: 
' the cosmological constant, the initial singularity problem, and the origin of the three spatial dimensions and time. 

, For the initial singularity problems, the tachyon condensation of winding strings may play a role when the radius of 
\^ ' universe shrinks smaller than the string scale, according to a rcccntproposal Here, we focus on the dimensionality 
\ problem in the context of the Brandenberger-Vafa (BV) scenario Q (see [1, In the original BV scenario 0, it is 
' assumed that all nine spatial dimensions start from the toroidally compactified radii near the string length and the 
'^j \ universe is filled with an ideal gas of fundamental matter, so-called string gas. It is also assumed that the string gas 
^ I is initially thermalized at the critical temperature Tjj, called the Hagedorn temperature In order to resolve the 
^ ' dimensionality problem, the string winding modes play a particularly important role. The winding strings prevent the 
, dimensions which they wrap from expanding, as shown in @, The annihilation of winding and anti- winding strings 
determines how many dimensions expand and how many dimensions stay at the string length. A simple counting 
argument suggests that this annihilation occurs mostly in the space-time dimensions of I? = 4, so the three spatial 

■ dimensions become large. Some studies have already examined various aspects of the BV scenario and it has been 
H ; extended in a variety of ways [1, 0, i, i, E [HI, E H Q [H, [H, [13 • (See also recent reviews 0, [H HI ) . 

The interaction (annihilation) rate F of the string winding modes, with the string coupling e*^, is roughly 



I. INTRODUCTION 



1001ni;e'''^. (1.1) 



(See [ill and Sec. |V] for more detailed discussion). Assuming that this interaction works efficiently at an early stage 
of the universe, decompactification proceeds. However, there arc three main assumptions on this scenario: adiabatic 
evolution, weak coupling and thermal equilibrium. A critical point of this scenario is the last assumption about 
thermal equilibrium [ll| (see also [l^l for another assessment). 
The thermal equilibrium condition is given by 

T> H, (1.2) 

where H is the Hubble expansion rate. Let us briefly estimate this condition based on a typical cosmological evolution 
at the Hagedorn temperature with "string matter" . The radius of the universe asymptotes to a constant value, whose 
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order is the string scale or several times of it. This means that all dimensions will be still small, which is inconsistent 
with our large four dimensions. However, if the radii (of three dimensions) asymptote to a sufficiently large value, 
then the universe could be matched to a radiation dominated phase. This type of evolution is only allowed when the 
initial dilaton satisfies a certain condition. From the condition indicated with Eq. (jl.ip . the interaction rate is found 
to be bounded from above as 



where the subscript denotes the value at an initial time (Sec. |V|. This shows that the thermal equilibrium condition 
is not satisfied at the initial phase of the universe, under the adiabatic condition H < 1. Therefore, no consistent 
treatment with well-defined thermodynamic functions can be carried out [2l] |. and all analysis based on such an 
approach cannot be trusted. 

This non-equilibrium aspect of string gas is related to cosmological expansion at the Hagedorn phase. Some features 
are necessary for resolving this difficulty, for instance: 

• The universe docs not evolve toward a constant radius at late times. 

• Initial evolution of the universe is modified. 

In the latter case, the universe may or may not evolve toward a constant radius, and the initial constraint for 
the dilation may be replaced by a weaker condition. There are several possibilities for achieving the above two 
evolutions. In the straightforward derivation of the interaction rate (|1.3p . we have assumed a simple dilaton-gravity 
system without taking into account any nontrivial coupling between the dilaton and "string matter" (gas) fields. Any 
nontrivial coupling modifies the dynamics of dilaton and cosmological expansion, and it may weaken the difficult 
requirement for the initial condition. For example, incorporating the NS-NS and R-R fields into the action p^. [l5|. 
or the effects of higher curvature corrections |16| will alter the evolution of the universe at an early stage. 

With these ideas in mind, we wish to pursue the possibility of resolving the thermal equilibrium issue. We adopt a 
simple modification of the scenario by taking into account the dilaton potential. This simple alternative will allow us 
to study the system rather extensively and give insight into other possibilities and approaches. We will see that the 
above two features are realized in the simple models in the present paper. 

The paper is organized as follows. In the next section, some general aspects of dilaton-gravity and string gas in 
the extreme Hagedorn regime of high-energy densities are briefly reviewed. We will take the dilaton potential to be 
of two types. The first case is discussed in Sec. HI, the second in Sec IV. For both cases, we analyze the dynamics of 
the Hagedorn regime with a single scale factor (Hubble radius) and with large and small radii. In Sec. V we discuss 
the thermal equilibrium of string gas for the two models. The final section is devoted to summary and discussion. 
We adopt the string scale a' = 1. 



Dilaton-gravity comes from the low-energy effective action of string theory. Ignoring contributions from the anti- 
symmetric two-form and including a potential V{(l)) for the dilaton, ^ the action of this system is described as 



where g is the determinant of the background metric g^i, and Cm is a Lagrangian of some matter. The coupling of 4> 
with gravity is the standard one arising in string theory. Assuming a spatially-homogeneous universe. 




(1.3) 



II. STRING GAS IN DILATON GRAVITY 



A. Basic equations 




(2.1) 



9 



(2.2) 



^ We give the dilaton potential in the string frame, while the corresponding potential in the Einstein frame is given in Appendix |A] Note 
that the initial data in the Einstein frame also differ from those in the string frame. 
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we can reduce the action to 

9 

S = / dtV^i-9'°) [e-^iY. - + - ^(^- PV^)] , (2.3) 



" 4=1 

where we have introduced a shifted dilaton ip, 



i^ = 20-^A,, (2.4) 



to simphfy the obtained equation of motion. In the reduced action, we have introduced the (one loop) free energy F 
of a closed string as the Lagrangian of matter. This is only possible in an early universe in which the string gas is 
in thermal equilibrium at the temperature /3~^. Variation with respect to goo, A^, ijj yields the following equations of 
motion: 



1=1 

K-i^\ = \e'^P^ + \v'{ct^), (2.5) 

i=l 

where E = pe^^' is the total energy and Pi = piC^^' is the total pressure along the i-tli direction, respectively, 
obtained by multiplying the total spatial volume e^^' with the energy density p and pressure pi. These quantities 
are related to the free energy F by the basic thermodynamic relation: 

dF 

Employing Eqs. (j2.5p . the conservation of the total energy is 

9 
i 

or equivalently, 

9 

p + ^A,(p + K) = 0. (2.7) 

i 

The basic equations become simple forms in terms of the shifted dilaton. In some cases, however, the equations 
written by the original dilaton are convenient. 

Here, if we separate the spatial dimensions into d-spatial large dimensions and (9 — (i)-spatial small dimensions, 
which are denoted as 

i? = e^, r = e'', (2.8) 
respectively, Eqs. (|2.5p in terms of the original dilaton take the following forms: 

- dfi^ - (9 - d)v^ + (20 - dfi - (9 - d)vf = e^-^p - V{<j>) , 

- (20 - d/i - (9 - d)i>)fi = ^e^^pd + \v'i^) , 

- (20 - d/i - (9 - d)i')i> = ^e^^p o-d + ^^'(0) , 

■■ -2 d(d-l).2 (9-d)(8-d).2 d(9 - d) . . 

4>-4>^+ ^ V ' + ^ J + ^ ^ > ^ 

= \e'^{dpd + (9 - d)p9_d) + \v{4>) + r(0) , (2.9) 
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where 

dF 



dv, 

in terms of the free energy F. 



P<^ = -T^e-''^-^9-'')^ (V* = l,...,d) (2.10) 
P9-d - -^e-''^-^'-''^'' , (V* = d + 1, . . . , 9) 



B. Hagedorn regime 

In the original BV scenario the very early universe is compactified over all nine spatial dimensions with radii 
R ^ 1 filled with string gas in thermal equilibrium around the Hagedorn temperature Q. In this regime, the 
usual thermodynamic equivalence between the canonical and microcanonical ensembles breaks down and the latter, 
more fundamental ensembles, must be used. Therefore the microcanonical density of states 0(£') needs to be derived 
by analyzing the singularities of the one- loop string partition function in the complex /3 plane, where /3 = 1/T is the 
inverse temperature. This analysis was carried out in [lol. [2^. Following these discussions, the partition function has 
poles depending on the radii of the universe, and accordingly the equation of state changes. 

For small radii, the leading order expression for the temperature and the pressure which come from the density of 
states is derived in [l3|. For large enough E, as a result, the temperature remains close to the Hagedorn temperature 
Th, and the pressure is vanishingly small, 

Pg-d-O. (2.11) 

Therefore, the string gas in the universe with small radii at the Hagedorn temperature can be treated as a prcssureless 
fluid, as expected by T-duality. 

As the universe grows with larger radii, the above mentioned poles move on the /3 plane. This causes the density 
of states J7 to change and yields a different temperature and pressure from those of the former "small radius" regime 
[2^. There will exist a critical radius R. When the radius is below the critical radius, the universe is still described by 
a pressureless state at the Hagedorn temperature, while above the value, the temperature and the equation of state 
are 

1 _ I3hE - 9 

where d denotes the expanding spatial dimensions. Note that these quantities depend on the total energy E and that 
it is necessary to solve the equation of motions for E as a function of time numerically. It was shown that the radius 
rapidly expands like an acceleration expansion while the dilaton continues its monotonic decrease [ll| . 

The above small and large radius phases are the basic states of string gas near the Hagedorn temperature, and the 
BV mechanism explains how the particular spatial dimensions enter into the second phase. These discussions are, of 
course, based on the assumption of the thermal equilibrium of string gas and that is the point of our discussion in 
this paper. 

When the energy density decreases and the temperature falls much below the Hagedorn temperature, the string 
gas will behave as radiation, with P = yielding the radiation dominated evolution, R oc t^/^'^+i) in [13, [HI- In 
what follows, we discuss the evolution of spacetime and dilaton and, based on the analysis, we study the thermal 
equilibrium of string gas at the initial Hagedorn regime, where the string gas behaves as pressureless dust. 



III. TIME EVOLUTION OF UNIVERSE AND DILATON: MODEL I 



We will take the potential of the dilaton to be of two types, each with a different nature. The first example is 
a simple exponential potential (Model I), while the second is a double- well potential, discussed in the next section 
(Model H). 

First, we shall discuss the exponential potential case. In order to describe the typical behavior of the system, we 
analyze the simple situation, in which all radii are the same, = v. The more generic situation in which small spatial 
dimensions (r = e") and large spatial dimensions {R = e^) are separated will be discussed later in this section. 
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A. Asymptotic solution with exponential potential 



As a toy model, we discuss a dilaton field with a runaway potential. (See [24| for early works in other context), 

V{(j)) =b e^"''' , (3.1) 

where a and b are constant parameters. For simplicity, we assume 0<O with a > hereafter to achieve the weak 
coupling regime e"^ < 1. Additionally, we take |6| = 1 since the basic behavior does not change according to the 
magnitude of b. The effective potential in the Einstein frame will allow us to intuitively understand the behavior of 
the dilaton. As discussed in Appendix El the effective potential of the dilaton in the Einstein frame W{<f>) is described 
as 



W = -e^V, (3.2) 

and hence the case with & < is the usual exponential potential in the Einstein frame. We will also analyze this case. 
For the following two limits, analytic asymptotic solutions can be found. Let us begin with discussing these solutions 
before we present numerical solutions. 



1/ — > limit 



This case is reduced to the standard analysis without the potential term. The equations (|2.5p with zero pressure 
give 



(3.3) 



(72/i^ ~ 36fi(t) + 4(/)2 + V) , 



and then the V approximation yields 



i', 



^1 • 



(3.4) 



dt- 



E 
2"' 



Since the total energy is constant. E 
and ^ are easily obtained from Eq 



Eo = 
as 



const, obtained from Eq. (|2.6p . the analytic solutions for 7/'(= 2(f> — 9/i) 



fi = fi^ + -\n 



{Eo{t - U) + 2B - 



-&A){B + iA) 



(3.5) 



{Eo{t - U) + 2B + 6A){B - 3A) 



where is the time when the condition V^(0) ~ is satisfied. The integration constants are given by A = /t,e' 
B = —ip^e~''^' , where and -0, are the field values at the time t ^ t^. 



-1/.. 



The solutions behave asymptotically in three ways (sec Fig. 1(a)). The first type (i) is that in which the radius 
grows as limt^o In IV^I ^-^^ diverges — > oo) at Eo{t — t*) = 2\B + 3A\. At the same time, the dilaton also grows 
and diverges ^. These behaviors appear if the conditions B + 3A < and \B + 3A\ < \B — 3A\ are satisfied, yielding 
< G/i* < (/)*. The second type of evolution (ii) is that in which the radius contracts as limf^o In |^| and diverges 



^ At the strong-coupling regime after the growth of the dilaton, wc cannot predict what will happen. A possibility of a static universe 
in the strong-coupling regime is recently discussed in llTll . 
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1: Evolution of asymptotic solutions p.Sp with respect to the initial data (fi,,<j),). Other initial data are fixed, as 
<j)t = with t, = 0. Fig. 1(a) shows the asymptotic solutions with V = 0. The forbidden region E < is equivalent to 



FIG. 

/X* = 

^ 101 ^ 6|/i|. The growing and contracting solutions correspond to < /i, < (/),/6 and /i« < </)*/3 with /i* < 0, respectively. 
For the other choice of initial data {(it,4>*), the radius converges to a constant value /Xoo of Eq. (|3.6fl . Fig. 1(b) shows the 
asymptotic solutions for the exponential potential V = he^"''^ with a = 2 and 6=1. The forbidden region E < is equivalent 
to (9/i — ^9/1^ — l)/2 ^ <^ ^ (9/i -f ^9/1^ — l)/2. The rapidly growing and convergent solutions correspond to 0c < 4>* and 
0, < 0c respectively. The critical velocity 4>c is approximated by straight lines of Eq. (|B9[I (/i, ~ (;/>*/3 + 5/6) with /i, > 
and Eq. (fBTTj) (/i. ~ 2/90*) with /i. < (Appendix B). 



at i?o(i — i*) = 2|i? — 'dA\. Similarly, the dilaton contracts. These behaviors appear if the conditions B — 'iA < 
and \B — 3A\ < \B + 3A\ are satisfied, yielding < 0, with /i» < 0. The last type (iii) is that in which the radius 
converges to a constant value, which is achieved for < B + 3A, B — 3A. 

The first two types of evolution wiU be seen in the double- well potential case (Sec. IIV[) . In the present model, only 
the last behaviors are important and we comment on its asymptotic radius. The asymptotic evolution at t ^ oo of 
the spacetime is 



B + 3A 



B-3A 



1/3 



(3.6) 



and it converges to a constant radius Q, while the (shifted) dilaton rolls monotonically to the weak coupling. Without 
fine-tuning {B ~ SA), the asymptotic value Roo is not very large, and the radius remains small. We call this solution 
the convergent solution. 

2. V,V' :$>(j)ij, 0^, fi'^ limit 
In the potential dominated case, the equations p.3p are reduced to 



0=iy(0) + y'(0)-6^e2-^ 



(3.7) 



As for & = T if the initial condition satisfies 



4a 

the analytic solutions for these equations are given by 



(3.8) 



8- 



1 



■ sm 



a I zGi V a 

M = ^^VT*^ + Di{t- t^) + L>2- 
8a + 1 



aCi{t-t, + C2) 



(3.9) 
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Parameter 




^ 'c — \ — '■ 

Type of solution 




> U 


a > 1 


convergent solution (00 < 0c) 








rapidly growing solution (0o ^ 0c) 


V = he^'"!' 




a < 1 


rapidly growing solution 




f)< 


a > 1 


convergent solution (0o < 0c) 








contracting solution (0o > 0c) 








convergent solution (iio > and 0o ^ 0) 




A > 
A < 


?? > 1 


rapidly growing solution (0o < 0) 

stabilized dilaton with ^ — ry. [Ha > and 00 < 0) 



TABLE I: Summary of basic behaviors 



Here Ci, C2, D\ and D2 are the integration constants given by 

C2 



a + 1 
1^' 



1 



aCi 

= A* - 

D2 = - 



■ cot 



2a 



2a 



1 



1 



(3.10) 



/i* and 0* are the field values at the time t = vifhen the condition V, V' 3> 0/i, 0^, /i^ becomes a good approximation. 
Both the dilaton and radius diverge at 



{t-t,)^—-C2. 

We call the solution (|3.9p the rapidly growing solution. Note that for 
Eqs. ([XT^ with 5=1.3 

For 6 = — 1, the solutions become 



< 



(3.11) 

b1 the solutions are given by 







1 


ln< 




i + - si 




a 



aCi{t-t,+C2) 



4a 



(3.12) 



The dilaton asymptotes to — Cit, toward the weak coupling, and at the same time the radii inflate if (8a + 
l)/t,/[2a(0, +Ci)] > 1; otherwise the radii contract exponentially. The contracting evolution appears in the numerical 
simulation since the condition for the rapidly contracting (8a + l)/i*/[2a(0* + Ci)]<l can be initially satisfied for 
|0o| ~ 0(1) and /io<0(l). 



B. Numerical analysis 

Now we are ready to present numerical solutions of the basic equations p.3p . We first discuss the case in which all 
spatial dimensions have the same radii, Xi = fi. The other case with Xi fj, will be discussed in the next subsection. 

In order to see the typical evolutions of dilaton and spacetime, we change the initial velocity 0o for various values 
of a, with respect to the fixed initial data of /io, 0o and fiQ (Appendix IB|) . We are interested in the situation where 



However, this solution does not appear in the numerical simulation in the next subsection because the dilaton goes to —Cit, violating 
the V-dominated condition. This case is, as a result, reduced to V — > limit in Sec. IIII A II 
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FIG. 2: Plot of typical convergent solutions of the dilaton <j} (blue line) and radius R = (red line) for /lo ~ (f'o = 0, fio ~ 
0.2, (j)o — —3 and a = 2,6 = 1. The analytic solutions (green and pink lines) of Eq. p.Sp are also plotted as a reference. Here 
the integration constants, ^, , 0*, /i*, are chosen appropriately. The asymptotic evolutions are well approximated by the 
analytic solutions. The same results can be obtained for the case <f)o > with /to ^ 0{1). 




i;- |i (analytic) 

I \. 10(1) 



1 (|) (analytic) 




all initial radii are near the string scale, R 1^ and where they expand (or contract) with time. Thus, it is natural to 
consider fiQ ^ and |/io| < 1, where the latter condition comes from the adiabaticity of the cosmological expansion. 
In our discussion, let us follow Table HI 



1. b > case 



Initially expanding case (Hq > 0) 

We begin with studying the initially expanding case, Hq = /ip > 0. For the initial condition of the dilaton, we 
consider (l>o<0 and (j)Q < 0, which satisfy the weak coupling condition. The case of 00 > will be discussed soon 
hereafter. 

a. Initial data of (/jq < 0. We set these initial conditions for various values of a and numerically solve Eqs. (|3.3p . 
For a > 1 wc find that there is a critical value of 0o, denoted (j)c- The critical velocity c/jc is approximated by a 



straight line in the (/io, ^o) plane, and the approximated line will be obtained by Eq. (|B9p in Appendix IB] (Fig. 1(b) ), 
where we will find (jjc < for Ij,q<0{1). For 0o < (f'Cj the radius goes asymptotically to a constant value, and the 

dilaton monotonically decreases, while both radius and dilaton diverge at late times for (f>o > 4>c- These two typical 
evolutions may be naturally related to the analytic solutions in the previous subsection. We plot these results in Figs. 
[5] and [31 comparing them with the analytic solutions. Fig. [2] shows a typical result for which corresponds 

to the convergent solution (|3.5p . and Fig. [3] shows the opposite case, which is related to the divergent solution (|3.9p . 
with appropriate choice of integration constants, {^*, 0*, /i*, 0*}. 

For any value of 00 < 0c j the behaviors of numerical solutions are the same as in Fig. [2l Similarly, for any value 
of 00 > 0c the behaviors are the same as Fig. [31 For 0o ^ 4>c, the radius seems to approach a constant value at 
first, but jl changes its sign from minus to positive at late times, and then /i eventually diverges. All these asymptotic 
evolutions are well approximated by the analytic solutions 

For < a < 1, the critical velocity and the convergent solution disappear, and all numerical solutions become the 
rapidly growing solutions, as in Fig. [31 



* It is also interesting to see if the numerical solutions satisfy the condition for the rapidly growing solution Eq. 1 13.81 1 at the initial time 
[t = 0). The condition Eq. 1 13. 81 is initially satisfied for /jo>0.3 with a = 2, while for fio<0.2 the numerical solutions initially violate 
the condition Eq. 1 13.81 1. but they immediately become the rapidly growing ones. 
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b. Initial data of (j)o > 0. At first sight, solutions with 0o > seem to easily violate the weak coupling condition. 
However, there are particular sets of parameters (/io,0o) that preserve the weak coupling condition. The same 
condition described above is also applicable in this case: for any value of < 0C: the numerical solutions are 
asymptotically the convergent solutions, as in Fig. [2l while for (j)Q > (j}c, the numerical solutions are asymptotically 
the rapidly growing ones as in Fig. [31 For 0o > 0, the convergent solution is allowed; otherwise, the weak coupling 
condition is immediately violated. This is realized by a positive critical velocity larger than (j)Q {(pc > 4>o > 0), yielding 
large initial velocity of radius Ao-^l (Fig- HI)- However, this case violates the adiabatic condition /i>0(l). 

Initially contracting case (Hq < 0) 

Next, we discuss the initially contracting case, Hq = /io < 0. In this case, we also find that the asymptotic behavior 
is the same as Hq > 0, even though the radius initially contracts. The critical velocity cjjc is also approximated by a 
straight line (see Eq. (|Blip ). The slope is different from that of Hq > 0, and we find (pc is negative in the present case 



(Fig. 1(b)). Hence for (f>Q < 0, any initial velocity of 0o < 0c yields the convergent solution, and the opposite case 
yields the rapidly growing one. On the other hand, for (j)Q > 0, all numerical solutions become the rapidly growing 
ones since 0c < 0. 

In summary, we find that for a > 1, there are mainly two asymptotic behaviors, i.e., either convergent or rapidly 
growing, depending on the initial dilaton velocity c/jq with respect to the other fixed initial values. On the other 
hand, for a < 1, we do not find evolution like the convergent solution. This feature will be understood clearly in a 
phase-space analysis in Appendix [Bl 



2. b <0 case 



Initially expanding case (Hq > 0) 

a. Initial data of (pQ < First, we consider the initially expanding case with 00^0 and 0o < 0. For a > 1, we 

find that there is a critical value 0c j similar to the case where 6 > 0. The asymptotic behaviors are, however, different 
from the previous 6 > case. For <f>Q < c/jc, the radius asymptotes to a constant value and the dilaton monotonically 
decreases, while both radius and dilaton contract exponentially for 0o > 0c ^■ 

These two typical evolutions are related to the analytic solutions of Eqs. (|3.5p and p.l2p in Sec. IIII Al Because 
0c is positive in the initially expanding case, the solutions become the convergent for all initial conditions 0o < 0. 



The exponentially contracting phase is not interesting from the point of view of decompactification; thus we will mainly pay attention 
to the convergent solution. 
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Notice that the forbidden region, E<0 for the b = —1 case, which is equivalent to the condition (9/io — 
a/1 + 9/j,g)/2<(/)o<(9/io + \/l + 9/ig)/2, exists near the critical velocity. So the critical velocity (pc is the same 
as either above boundary value of this forbidden region. The critical velocity (j)c is also approximated by a straight 
line in the {fiQ, 4>o) plane as in the case where b > 0. The slope is different from 6 > and the above forbidden band 
E<0 is near the line. 

h. Initial data of (jjo > The same condition as above is also applicable in this case: for any value of (t>o < 4>c^ 
the behaviors of numerical solutions are asymptotically the convergent solutions, while for any value of 0o > 0c j 
the numerical solutions are asymptotically the rapidly contracting ones. Although 0o > seems to violate the weak 
coupling, the weak coupling condition is held asymptotically since the dilaton goes to the weak coupling region for 
both the convergent solution and the contracting one. 

Initially contracting case (Hq < 0) 

Next we discuss the initially contracting universe. In this case, we also find that the asymptotic behavior is the 
same as Hq > 0, even though the radius initially contracts. The critical velocity (jjc in the present case is also 
approximated by a straight line, and we find (f>c is negative. Hence for (f>o < 0, any lower value 0o < (f'c yields the 
convergent solution and the opposite case yields the rapidly contracting one. On the other hand, for <^o > 0- 
numerical solutions become the rapidly contracting ones because (j)c < 0- Notice that the critical velocity 4>c' is also 
the same as either boundary value of the forbidden region of E<0, similar to the case where Hq > 0. 

In summary, we find that for a > 1 (and b < 0), there are mainly two asymptotic behaviors, i.e., the convergent or 
rapidly contracting, depending on the initial dilaton velocity (j>Q with respect to the other fixed initial values ^. In the 
initially expanding case, for 0o^O, all numerical solutions become convergent solutions. The behavior will be easily 
understood in the effective potential picture in the next subsection. 

3. Effective potential picture 

As mentioned above, the effective potential of the dilaton in the Einstein frame is = — 6e^'"^ with a redefined 
constant d = a + 1/4 obtained from Eq. (|3.2p . Adopting the Einstein frame, we can explain intuitively why the 
evolutions depend on an initial velocity (j)Q. 

b > case 

The effective potential has a valley in the strong coupling regime (cj) > 0) and asymptotes to a flat value, W ^ 0, 
in the weak coupling regime {(j) < 0). In this picture, when the dilaton starts from the origin (/)o^0i with a small 
initial velocity, it rapidly falls down a positive side of the potential. On the other hand, if it starts with a sufficiently 
large velocity to avoid falling, it can climb up the valley of potential toward W = 0. Therefore, V becomes dominant 
for small initial velocities, while for large initial velocities, V goes to zero. These behaviors are consistent with the 
above-mentioned two regimes in the analytic method. 

6 < case 

This case is equivalent to the usual picture of exponential potential in terms ofW"^, where the dilaton rolls toward 
W ~ in the weak coupling. In this picture, when the dilaton starts from the origin 4>o<0, it tends to go to F = 0, 
asymptoting to (p — oo. For the initial condition, 0o < or < (/)o<0(l), the solutions become V ^ 0, while for 

(/> > 0(1), fi begins to rapidly contract since the condition (8a + l)/i*/[2a(t/)* + Ci)]<l is satisfied before V goes to 
zero, as shown in the previous discussion. 



® The classification of the solution for a < 1 is more complicated. For the initially expanding case with 1/2 < a < 1, in addition to 
the above two evolutions, there is the rapidly growing solution for the very large velocity (/)o 2> 1. As for a < 1/2 the contracting 
solution disappears. Hence the asymptotic solutions are of two types; the convergent and rapidly expanding ones. However, the rapidly 
expanding solution in the present case initially violates the weak coupling condition because <po > 0. On the other hand, for the initially 
contracting case with a < 1, there exists convergent and contracting solutions similar to those where a > 1. This a-depcndence will be 
clarified in Appendix [B] by employing phase-space analysis. 

As a natural potential for the dilaton, W oc e'^* is interesting to study [25l. l26l|. It physically corresponds to the case a 3> 1 with b = —1 
in our setup. However, the basic behaviors do not change and are similar to the case where a > 1 described above. In fact, we have 
confirmed that the basic behaviors of the dilaton with W = e'^'* are similar to the case where a ^ 1. 
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FIG. 4: Typical convergent solutions of radii of large spatial d = 3 dimensions {R — e'') and small 6 dimensions (r = e'^). The 
initial conditions are no = uq = <j)o — 0, fio — 0.8, !>o = 0.1, (po — —3 and a — b = 1. 
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FIG. 5: Typical rapidly growing solutions of i? = e'' and r = e"^ for the same initial conditions and parameters as in Fig. |4] 
except for 00 ~ —1. 



C. Evolution of large and small radii 

So far, we have assumed all spatial directions have the same radius. As expected, the basic behaviors observed 
under such simplification hold even when there is asymmetry between several spatial radii. We briefly discuss such 
cases, representing a typical example of evolution with different spatial radii. 

b > case 

In Figs, m and O we plot typical two types of time evolution of radii, as a solution of (|2.9p for d ~ 3 with the 
potential (|3.ip . From these figures, it is understood that each radius evolves in the same way: both of them are 
either convergent solutions or rapidly growing ones. Wc could not find numerical solutions in which the radii show 
different types of behavior in late times. The reason why such solutions arc not found is as follows. Suppose the 
potential V dominates in the equation of motion for /x, then v becomes the convergent solution asymptotically. This 
assumption makes /i rapidly grow. In this case, the derivatives of v are asymptotically vanishing, and by substituting 
this condition into the second of Eqs. (|2.9p . we obtain the equation V' ~ 0. This is inconsistent with our assumption 
that V is dominant. Therefore the rapidly growing solution of fi is inconsistent with the convergent solution of and 
vice versa. 
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We conclude that for a > 1 the behavior of each radius is determined by the initial velocity (f>o once fiQ and i^o are 
fixed. Both radii evolve as the rapidly growing solutions for a large dilaton velocity while they are the convergent 
solutions for a small velocity. 

& < case 

Both fi and z/ behave in the same way as they do in the case where b > 0: both radii are convergent solutions or 
rapidly contracting ones since the same argument for & > is applied in the present case, except for the different 
evolution during the ^-dominated regime. For a > 1, 0o < and < liQ,i>Q<0{l), both radii always become 
convergent solutions, and the typical evolutions are the same as in Fig. [H 



IV. TIME EVOLUTION OF UNIVERSE AND DILATON: MODEL II 



A. Double-well potential and asymptotic solutions 



As a second model, we consider the double-well potential 

A 



(4.1) 



where A and 77 are a coupling constant and a vacuum expectation value (VEV) of 0, respectively. Similar to the 
previous model, A is taken to be positive or negative. The later case corresponds to an ordinary concave effective 
potential picture W ~ —V (Appendix E| . For simplicity, we assume until jjlV CI that all radii are the same, n = v. 

For A < ®, the dilaton stays at the minimum (VEV) of the double- well potential (pv = ±??. Expanding cj) around 
the VEV, 



(4.2) 



we find that Eqs. (|3.3p are approximated as 



/i 9/i2 w , 

77^ 4 



(4.3) 



These equations have the following solutions: 



£1 
9 



(4.4) 



?] \ 9 

where Ju{z), N^{z) are the Bessel functions with the amplitude Ci and C2. Di and D2 are given by 

1 



Di = — 



D2 = (U* 



ln9|/t* 



(4.5) 



Note that when we derive the above solutions, we have used the fact that the last term in (|4.3p decays as e^^^^^"^^' ex 
1/t, so that it can be omitted at late times. The amplitude of 54> becomes cx 1/Vt for A < 0, which is equivalent to a 
decaying solution. We will call this solution the stabilized dilaton solution. 



For the case where A > 0, we find solutions like Eqs. 114.41 1. which represent the dilaton oscillating around <f) = 0. 
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B. Numerical analysis 



Now we are ready to present the numerical solutions of the system. We take the initial condition in such a way 
that /io(~ 0), (/)o(<0) and /ig ^.re fixed and 4>o is changed for various values of A and ry, within the adiabatic condition. 
The evolutions depend on the sign of A, and so we begin with the case A > 0, following Table HI 



1. A > case 

We find that the potential is negligible and hence does not work for the large initial velocities > 3/io and 
t>o > 6/io with the initially expanding and contracting case, which are equivalent to the regions labeled as "growing" 



and "contracting" in Fig. 1(a) For these initial conditions, the numerical solutions are the same as the analytic 
solutions (j3.5p with = to, obtained by y ^ approximation. As discussed in §111. A, the asymptotic behaviors are 
two types: both the radius and dilaton either grow or contract. However, from the point of view of decompactification, 
the contracting solutions are not interesting. Moreover, the growing solutions violate the weak coupling condition 
under 0o ^ due to the positive velocity, 4)o{> 6/io) > 0. Therefore, we discuss the other initial conditions, 0o < 3/io 
(or 00 < 6/io), for which the potential works. 

In this case, there exists a critical value of (/jq, denoted 0c, which determines the evolutions as follows. The radius 
begins to contract at late times as the dilaton rolls to the negative infinity for 0o < 0c ? while the radius diverges at 
late times as the dilaton rolls to the positive infinity for 0o > 0c- This behavior is easily understood by employing 
effective potential picture. The effective potential of the dilaton is given hy W ^ —V. In this picture, for A > 0, there 
exist two valleys in both sides of the potential for |0| oo. When the dilaton starts from the origin, it asymptotically 
goes to the positive (negative) infinity, |0| — > oo, if the initial velocity is 0o > 0c (0o < 0c)- On the other hand, 
taking 0, 1/ ^ 1 limit yields V ~ A/4(0/r7)^. Combining this with the fact that the total energy is constant from 
(|2.6p . E K, ^i^-'^'^y ~ const, we obtain e^^ oc e^'^0~'*. This equation means that the dilaton going to the positive 
(negative) infinity makes the radius go to the positive (negative) infinity (|/i| oo). /i oo corresponds to the 
rapidly growing solution, and // — > — oo corresponds to the contracting solution. 

Even though the critical velocity 0c depends on both parameters A and r/, it depends more sensitively on 77 than 
A. The value of 0c becomes smaller as A or 77 grows. At 77 1, 0c is positive ~ 0(1). Then for 77<1 with 0o < 0, in 

which the weak condition will be satisfied, the radius always begins to contract since the condition 0o < 0c is satisfied. 
Since we are mainly interested in the expanding phase, let us consider 7/ > 1. For a large value of 77(^ 1), e.g., 77>5, 

both the radius and dilaton expand and diverge at some late times for any value of 0o < since 0c (< ~O(30)) is 
sufficiently small to satisfy the condition for the rapidly growing solutions, 0o > 0c- It means the dilaton goes to the 
positive infinity for a very wide range of parameters. We will revisit this behavior soon later (Sec. IIVB 3[) . 

In summary, we find that there arc mainly two asymptotic behaviors, i.e., contracting or rapidly growing, depending 
on the values of initial dilaton velocity 0o. In the case 0o < 0, in which the weak coupling is satisfied, all numerical 
solutions become the rapidly growing ones for both the initially expanding and contracting cases as long as 77 ^ 1. 



2. A < case 



Similar to A > 0, the potential does not work for the large initial velocity of the dilaton, 0o > 3/io and 0o > 6/io. 
As discussed before, these behaviors are not suitable for our expected situation, and so we examine the other initial 
conditions where the potential works. In such case, the new evolution comes out: the radius grows monotonically, 
while the dilaton asymptotes to the VEVs (0y = ±7/), settling down to the potential minimum (Fig. [6]). The critical 
velocity 0c determines the VEV (0^ = —77 or 77), where the dilaton is stabilized in the end. For the initial velocity of 
00 < 0c J the radius monotonically grows and the dilaton asymptotically goes to the negative VEV —77, while for 
00 > 0Cj the dilaton asymptotically goes to the positive VEV 77. The former case corresponds to a stabilizing 
dilaton at the weak coupling. In Fig. [6l we plot the numerical results and the analytic solution Eq. (|4.4|) with 
appropriate integration constants, /i*,/i,t. Note that the forbidden region E<S) exists near the critical velocity 0c in 
the present case (Fig. [1} and then that the allowed region of 0o which makes the dilaton go to the positive VEV 
(0v/ = 77) is narrow since the larger value of 0o yields the "growing solution" indicated in the F ~ case of Fig. [T] 
Therefore, the dilaton goes to the negative VEV for a very wide range of parameters. In the initially expanding case, 
for any value of 0o < 0, the radius monotonically grows and the dilaton asymptotes to the negative VEV because 0c 
is positive for Hq > 0. The result of these behavior docs not depend on 77. 
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FIG. 6: Typical evolutions of dilaton and radius for the stabilizing solutions with the same initial conditions as in Fig. |3]except 
for A = — 1,»7 = 1. The analytic solutions (|4.4|l are also plotted as a reference. The dilaton goes to the VEV {(f>v = —rj) as 
oscillating. The radius expands monotonically. 



3. Effective potential picture 

At first sight, the time evolution of the dilaton under the double- well potential may seem to contradict our intuition. 
As mentioned above, however, the effective potential of the dilaton in the Einstein frame is ~ —X{{4>/rff' — l)^e'^/^ 
from Eq. (|3.2p . Therefore, for A > 0, there exist two valleys on both sides of the potential \4>\ oo, and the slope of 
the valley existing in the strong coupling regime (0 > 0) is steeper than in the weak coupling regime by a factor e*^/^. 
Notice that the dilaton falling down to the positive side of the potential corresponds to the rapidly expanding radius, 
while the dilaton falling down to the negative side corresponds to the rapidly contracting radius. Hence when the 
dilaton starts from the origin, it is easier to make the dilaton roll into the positive side and fall down. Then the radius 
grows rapidly for a very wide range of initial velocities (j)Q. This observation is consistent with the above numerical 
results. 

As for A < 0, it is basically equivalent to the usual double- well potential picture in terms of W except for its 
different slope by a factor e"^/^, where the dilaton rolls toward the VEV. In this picture, when the dilaton starts from 
the origin with (j>o<Q and 0o<O, the dilaton always goes to the negative VEV, which is again consistent with the 
above numerical results. 



C. Evolution of large and small radii 

Similar to the exponential potential case, we finally discuss evolutions with different radii in the double- well potential 
case. 

For A > 0, as a result, each radius evolves in the same way: the contracting solution and rapidly growing one 
described as in Fig. [51 We cannot find a solution in which each radius evolves differently. The asymptotic behavior of 
the radius is determined by the behavior of the dilaton, i.e., rolling down to the positive side of the valley or rolling 
down to the negative side of the valley in the effective potential picture, which corresponds to the rapidly growing 
radius and the contracting one, respectively. Therefore, both radii evolve in the same manner since the potential acts 
on each radius in the same way, which is clearly seen in Eqs. (|2.9p . For (po < and |/io|, l^'ol^O(l), both radii always 
become the rapidly growing solutions as long as 77 3> 1. 

As for A < 0, similar to the above case of A > 0, both radii evolve in a same manner: the contracting, growing, 
and stabilized dilaton solutions. We plot the typical evolution of stabilized dilaton solution in Fig. [T] For 4>o < and 
< /xo, vo<0{l), both radii always become stabilized dilaton solutions, and such behavior does not depend on 77. 
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FIG. 7: Time evolution of radii large spatial d = 3 dimensions {R — e^) and small 6 dimensions (r — e") for the typical 
stabilized dilaton solutions. The initial conditions are the same as in Fig. [S] except for A = — 1 and r; = 0.5. 



V. THERMAL EQUILIBRIUM OF STRING GAS 

In order to resolve the dimensionality problem, the annihilation of winding strings play a critical role in the context 
of the BV scenario. The interaction rate F of annihilation, equivalently, the process where any winding/ anti- winding 
string pair annihilates to the momentum string pair with the coupling given by e"*, is roughly [ll[ 

r~ 1001ni;e'''^, (5.1) 

where 100 is a numerical factor of sum over both spins and momentum states and E is the total energy of string gas. 
The thermal equilibrium condition is given by 

V>H, (5.2) 

where H is the Hubble expansion rate. The assumption of the thermal equilibrium of string gas is necessary for the 
BV scenario to work. However, a typical cosmological evolution around the Hagedorn temperature is that the radii 
of the universe asymptote to a constant value i?oo as described by Eq. p.6p . As discussed in Sec. IIII Al Roo does not 
become large without fine-tuning. It means that all dimensions are still small and it is inconsistent with our large 
four dimensions. Nevertheless, if the radii (of three dimensions) asymptote to a sufficiently large value, the universe 
could be matched to a radiation dominated phase, and only such a case is a viable scenario. Therefore, we require 
that the asymptotic radius will be larger than a critical radius R, i.e., i?oo > R- Here i? is a characteristic scale which 
divides a "small" radius from a "large" one (Sec. Ill Bp . Then, from (|2.5p and (|3.6p . we obtain a constraint equation 
on the initial value of dilaton as 



Ei3 



{Re-t'^y + 1 



(5.3) 



where the subscript implies the value at an initial time. Substituting this initial constraint into ()5.ip . the interaction 
rate at an initial state is roughly estimated as 



r < lOOlns(^) O(lO-i) < (o(10-i)Fo^)i?o, (5.4) 

where we have taken i? 3, /ig ~ and sufficiently high energy, E ^ 100. This shows that the thermal equilibrium 
condition is not satisfied at an initial condition as long as the adiabatic condition Hq ^ 1 holds. In what follows, we 
test the thermal equilibrium condition (|5.2p for our models. 
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FIG. 8: Plot of T — H for the model of exponential potential {b = 1). The initial conditions are /io — 0,(t>o = —0.1, /io — 0.3 
and <f>o — ~1 (the rapidly growing solutions as red and green lines), or (po = —3 (the convergent solutions as blue and pink 
lines), r — H > means thermal equilibrium. For the rapidly growing solutions, the weak coupling condition (e'* < 1) breaks 
down at t > 3 for a = 1 and t > 4 for a = 0.1. 



A. Exponential potential case 

1. b > case 

The exponential potential with b > and a > 1 allows two types of trajectory for its solutions, i.e., the rapidly 
growing and convergent ones, described by the analytic solutions (j3.5p and p.9p . respectively. Firstly, the convergent 
solution makes the radius grow toward the asymptotic value of Eq. (|3.6p , yielding the constraint of Eq. (j5.3p . However, 
once the potential term is included, the constraint (|5.3p does not directly restrict the initial value of dilaton at t = Iq, 
but it provides a constraint at a late time, t = t^,. Therefore, it could relax the constraint on the initial dilaton value, 
resulting in better realization of the thermal equilibrium. Figure [5] plots the thermal equilibrium condition with time: 
r — i7 > is equivalent to the thermal equilibrium condition. From this figure, we find the thermal equilibrium is 
realized until t ^ 1, which is equivalent to the time scale on which the potential term works, t<t^. This duration 
does not depend sensitively on a. The final scale of the radius is, however, not so large, and it may be insufficient to 
continue the decompactification process. Nevertheless, a noteworthy point is that the thermal equilibrium is realized 
at the initial phase, contrary to the naive scenario. 

Secondly, the rapidly growing solution is the solution in which the radius grows unboundedly, and thus the initial 
constraint cannot be applied in this case. Contrary to the convergent case, the dilaton grows toward the strong 
coupling, and the weak coupling condition breaks down. We see that the thermal equilibrium is also realized as seen 
in Fig. [HI So the thermal equilibrium is at least realized until the time of the violating, even though the time of the 
breaking becomes longer as a decreases. The duration of thermal equilibrium becomes large at most by a factor of 
two. For a < 1, all solutions are the rapidly growing type, and, for example, the duration is about t ^ 5 for a ~ 0.01. 

2. b <0 case 

In this case, all solutions with a > 1 that satisfy the weak coupling condition become convergent solutions. Thermal 
equilibrium is the same as the convergent solutions in 6 > 0, and hence the initial equilibrium continues until t ^ 0(1). 
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FIG. 9: Plot of r — _ff for the model of double-well potential. The initial conditions are /io — 0,00 ~ — 0.1,/io ~ 0.3 and 
(/)() = —1. For the rapidly growing solutions, the weak coupling condition breaks down at t > 10 for A = 1 and t > 15 for 
A = 0.5. 



B. Double- well potential case 

1. A > case 

All solutions with 77 3> 1 that satisfy the weak coupling condition become rapidly growing solutions. In Fig. [5] 
we plot the thermal equilibrium condition. The basic behavior of thermal equilibrium is the same as in the growing 
solution for the exponential potential with 6 > 0, except for the time scale: the initial equilibrium is realized until 
t ~ 0(10) where the weak coupling condition violates. The thermal equilibrium is only marginally satisfied at almost 
all times until the violating. 



2. A < case 



In this case, all solutions which satisfy the weak coupling condition become the stabilized dilaton solutions. If we 
take rj to be small, e.g., 77 < 0(1), the thermal equilibrium continues unboundedly (Fig. This is because the solution 
stabilizes the dilaton at the VEV {(j> = —77) and then the interaction rate F oc e^'^ asymptotes to e"^** ^ 0(0. 1). At 
the same time, the Hubble expansion rate asymptotes to zero, H ex 1/t — > 0, and then we have (F — H)\t^oo > 0. 
As Fig. [5] shows, the thermal equilibrium continues unboundedly as long as ?7<1. This result is the best situation 
for the BV scenario in our models. As for 77 > 1, the thermal equilibrium is only marginally satisfied at late times, 
T-H^O. 



VI. SUMMARY AND DISCUSSION 



We have studied the thermal equilibrium of string gas in the Hagedorn regime where the universe is in high energy. 
Thermal equilibrium is the one of the important assumptions for the BV scenario in order to induce the dynamical 
decompactification of three large spatial dimensions. However, the initial thermal equilibrium condition of string 
gas is not realized in the original scenario based on dilaton-gravity. To resolve this difficulty, we have explored 
possibilities for avoiding the issue. As a first step to tackle this problem, we have studied a minimal modification of 
the original model, by introducing a potential term of the dilaton. This simple setup allows us to study the system 
rather extensively. However, this does not mean that stabilization of the dilaton or the effects of the potential term is 
a necessary ingredient. We wish to emphasize that effects of matter (e.g., flux or any kind of corrections, etc.) would 
not be negligible, and taking into account such effects, we could avoid the issue of thermal equilibrium in the early 
universe. We expect our simple setup will provide implications for such effects. 

We have taken the dilaton potential to comprise two simple potentials, i.e., the exponential potential and the 
double-well potential, and have analyzed both the dynamics of the system and the thermal equilibrium condition at 
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the initial stage of the universe. Even though we have mainly studied the evolution of the scale factor with same radii, 
there is no significant difference in the typical evolutions of different radii. Based on the solutions, we have examined 
whether they satisfy three basic assumptions, i.e., the adiabatic condition Ho<0{l), the weak coupling condition and 
the thermal equilibrium condition T > H . As a result, we find the following cases. 

Exponential potential : V — b e^^"^ 

(i) b > : the convergent and rapidly growing solution 

The convergent and the rapidly growing solutions for a > 1 are acceptable cases. The former case is that the radius 
converges to a constant value and the dilaton rolls monotonically to the weak coupling regime. In the latter case, both 
radius and dilaton asymptotically grow in short time. These different evolution is determined by whether the initial 
velocity of dilaton (f>Q is smaller or larger than the critical velocity (jjc- On the other hand, for a < 1, all numerical 
solutions are asymptotically the rapidly growing solutions. Both of two solutions satisfy the thermal equilibrium 
condition during some initial time. 

(a) b < : the convergent solution 

In this case, the convergent solutions are sensible solutions since the others are contracting evolutions. For any value 
of a, all numerical solutions satisfy the thermal equilibrium condition during the initial time. 

Double-well potential : V = |((^)^ - 1)^ 

(i) X > : the rapidly growing solution 

For 77 ^ 1, all numerical solutions are rapidly growing solutions, and they satisfy the thermal equilibrium condition 
during some initial time which is relatively longer than the above cases. 

(a) A < .' the stabilized dilaton 

In this case, the dilaton is stabilized and it does hold the weak coupling condition. For any value of r/, all numerical 
solutions are reduced to the stabilized dilaton solutions. If we choose 77<0(1), the thermal equilibrium continues 
unboundedly. 

From these results, we conclude that it is possible to realize the thermal equilibrium of string gas at the initial 
Hagedorn regime. At the end of the paper, we would like to ask if the evolutions can be matched to a late-time 
universe. For the convergent solutions, the universe will not enter into the large-radius phase R > R, because R does 
not become large without fine-tuning, as discussed above. Besides, the time scale of the thermal equilibrium will too 
short for the BV mechanism to work, compared with Hubble time. Similarly, for the rapidly growing solutions, the 
short time scale will be a problem in the case of the exponential potential. However, in the double-well potential 
case, the duration becomes relatively longer and the situation is better. Nevertheless, it remains as a problem that 
the dilaton asymptotes to the strong coupling, where we cannot predict what will happen. 

On the other hand, the stabilized dilaton solutions can make the radius grow as i? oc and it would be possible to 
match it to a late-time universe. Moreover, the thermal equilibrium continues unboundedly. This solution is the best 
case among all examples presented in our paper. It implies that the (quasi-) stabilized dilaton at an early stage of the 
universe improves the situation. Besides, this example implies an interesting possibility for constructing a model that 
resolves the stabilization and dimensionality problem at the same time. It remains an interesting question whether 
we can build a model which can be embedded in string theory resolving all problems described above. 
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APPENDIX A: EFFECTIVE POTENTIAL 

In this Appendix we will show the effective potential of the dilaton in the Einstein frame derived by conformal 
transformation. Quantities in the conformal metric will be denoted by a tilde. We will denote the conformal factor 
by £7^ which is a function of the D dimensional spacetime coordinates a;^. The conformally transformed metric is 
= ^'^Qfiu 1 and the determinant of the metric scales as \/—g = U^y'—g . 
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(Al) 



Consider the dilaton-gravity action in D dimensions, 

Under the conformal transformation, the action becomes [l^ 

[ d^xy^\e-^'^fl^-'^ \r + 4(V(A)2 + (£> - 2){D - + Q-^V + 4{D ~ l)n-\\/n){^(j)) . (A2) 



The last term comes from integrations by parts. If we set the conformal factor as ^ = e^"^ , the action in the 
Einstein frame is 



S = / d^x-vO^ 



(A3) 



In the case D = 10 the Lagrangian of the dilaton which is minimally coupled to the metric in this frame, is reduced 
to 



with the canonically normalized potential 



W{(j)) =-e^V{ct)). 



(A4) 



(A5) 



APPENDIX B: AUTONOMOUS PHASE PLANE 

In this Appendix, we describe the asymptotic behavior of the system with the exponential potential. We define the 
dimensionless phase-space variables [28l . |29| 



(Bl) 



2</>-9i7 J\V\ y^e'f 

x= — — — , y ^ „ „ , z = 



3H 



3H 



3H 



The Eqs. (|3.3p are reduced to a so-called plane autonomous system, which consists of three evolution equations and 
a constraint equation, 



-[I - + {1 - a)y^ ~ Saxy''] , 
' 3 

9 2^ [(2 - + 3a2/^ - 3a] , 
3 

--z{x + 3ay'^) , 



l = x'' +y 



2 z\ 



(B2) 



Here a prime denotes a derivative with respect to the number of e-foldings, ln(i?) = /i. As the system is symmetric 
under the reflection y — + —y and z —^ — z, we mainly consider only the upper half-planes, y > and z > 0, in the 
following discussion ^. We solve this system numerically to find the phase-trajectories. In Figs. [10] -El we show 
phase trajectories on the y-x and z-x planes for a — 1/2, 1, 2, 3. In this system, there arc three types of fixed 
point (critical point), defined by x' = 0,y' = 0, z' — 0: 



X = ±1, y = 0, z = , 



A± , ., 

2- a , 2V2a2 + a - 1 
B± : x^— — , y = ± z , z = U. 



3a 



3a 



^ y < corresponds to the contracting evolution (fi < 0), and the hnes of phase-trajectory are the same as y > 0, except the direction of 
the arrow (Fig. IIOI I. 
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FIG. 10: The phase trajectories for a = 2 in {x,y) and {x,z) planes. There are four fixed points in the upper half-planes. B+ 
is the attractor point for the rapidly growing solution. C+ divides the convergent solution from the rapidly growing one. The 
red lines show the trajectories approaching the rapidly growing solutions, and the black lines show the trajectories approaching 
the convergent solutions. The arrows represent the direction of time evolution. The color of the lines indicates the end state 
of these solutions. An enlarged figure showing the behavior of the black lines in {x, z) plane is depicted in Fig. 1111 The dotted 
lines represent the constraint 1 — x'^ + or 1 = x'^ — , which means the forbidden region in the phase space. 



C± : X 



3a 



z = ± 



3a- 2 



(B3) 



C± is vanishing for < a < 1, and B± is also vanishing for < a < 1/2. For a < 1, there arc four fixed points 

and (7+ in the regime y >Q and z > 0. For 1/2 < a < 1, there are three fixed points, A^, and _B+. For 
< a < 1/2, A+ and A_ are the fixed points. It can be understood from the global behaviors on the phase-space 
that A± may correspond to an unstable node (or saddle node). 

In the following we show that the critical point describes the rapidly growing solution of Eq. (|3.9p and C+ 
divides the convergent solution (|3.5[) from the rapidly growing one. In order to see the behavior of these solutions 
on the phase-space, we rewrite these solutions in terms of the phase-space variables {x,y,z). From Eqs. p.9p . the 
rapidly growing solution with b = ±1 at late times {fi — > oo) gives 



20-9/i 1 

X = — — cx - , 

3/i 3 

y = -— oc 1/ cos aC\{t — + C2) 
3/i 



sin 




1- 


1/0 


■^^ cos 







(B4) 



Therefore, for a > 1, x and y asymptote to finite values, while z goes to zero at the divergent point aCi{t—t^ + C2) ^ tt, 
as seen in Eq. (|3.1ip . The trajectories approach the attractor point which arc plotted by the red lines in Figs. [10]- 
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FIG. 12: The phase trajectories for a = 3 in {x, y) and (a;, z) planes. The four fixed points are the same as in Fig. 1101 



[T2l Similarly, all trajectories for 1/2 < a < 1 approach the attractor point (Fis. fT3|) . while for < a < 1/2 the 
critical point plays a role of corresponding to the rapidly growing solution (Fig. 

On the other hand, the behavior of the convergent solution at late times, when the (shifted) dilaton decreases to 
the negative infinity, ?/'(= 20 — 9fi) — > — oo, follows from Eq. (|3.5p : 



. = |^oc-f(.-t.)-B, 
3/i 

Therefore, x — *■ — oo and z ^ oo for a > 2, while y goes to a finite value at late times. In Figs. [TUlfT^ the trajectories 
that approach asymptotically the convergent solution are plotted by the black lines. 

From these figures for a > 1, we see two types of asymptotic trajectory: the trajectory approaching i?+ which 
represents the rapidly growing solution and the one approaching the convergent solution (|a;|, \z\ — > oo). These two 
types of trajectory are clearly divided by flows around the fixed point C-^.. For < a < 1, all trajectories approach 
the rapidly growing solution because the fixed point (7+ disappear as shown in Figs. 1131 and [Til 



22 




FIG. 13: The phase trajectories for a = 1 in {x,y) and {x,z) planes. j4+, A- and B+ are the fixed points. (7+ does not exist, 
contrary to the other cases, and hence there is no trajectory approaching convergent solutions. AH trajectories approach B+, 
and all these solutions are asymptotically rapidly growing solutions. 



y 




FIG. 14: The phase trajectories for a = 1/2 in {x,y) and {x,z) planes. A+, A- are the fixed points. B+ and C+ do not exist. 
A+ is the attractor point replaced by B+ corresponding to the rapidly growing solutions. All trajectories approach A+, and 
these solutions are asymptotically the rapidly growing solutions. 



From these figures, for a > 1, we can estimate the critical line which divides the rapidly growing solution from the 
convergent one in the phase-plane {x,y). The critical line is approximated by the straight line passing through yl_ 
and C+ 



yc,-^^ix + l). (y>0). (B6) 
2a + 1 



Therefore, the trajectories that asymptote to the convergent solution are given by 



1^(^ + 1), (B7) 

while the trajectories that asymptote to the rapidly growing solution are characterized by the opposite inequality 
sign. This condition can be rewritten in terms of the initial dilaton and spacetime variables. Assuming (j)Q ~ 0, the 
condition for the convergent solution is 

^o<3mo--^^, (B8) 
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which is appKcable for the expanding case fiQ > {y > 0). 

In our numerical analysis in Sec. IIIIBi we vary the initial velocity of the dilaton, fixing other initial conditions. 
From the above equation, the critical velocity for the initially expanding case that divides the late-time evolutions is 
approximately given by 

2a + 1 , , 

~ 3Ao - . . J ■ (B9) 

Any velocity lower than the critical velocity yields the convergent solution, while any velocity beyond the critical 
one yields the rapidly growing one. The analytic estimation (jB9[) is in good agreement with the critical velocity 0c 
obtained from the numerical analysis. For example, we find 4>c I'Pc ^ 0-88 for /io = 0.2 and (pc I'Pc ^ 1-3 for /to = 0.5, 
based on our numerical simulation for a = 2. 

For the contracting case y <0, the trajectories asymptoting to the convergent solutions are obtained from Fig. [TOl 

x>Q. (y <0). (BIO) 
The condition required for the convergent solution is 2(j)Q < 9/io, and then the critical velocity ipc for /io < is 

- 9/2/io. (Bll) 
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